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Abstract. In the present paper we introduce a notion of G— decompositions of matrices. 
Main result of the paper is that a symmetric matrix A m has a G— decomposition in the class of 
stochastic (resp. substochastic) matrices if and only if A m belongs to the set U m (resp. U m ). 
To prove the main result, we study extremal points and geometrical structures of the sets U m , 
' U m . Note that such kind of investigations enables to study Birkhoff's problem for quadratic 

G— doubly stochastic operators. 
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4^ 1. Introduction 

Let us recall that a matrix A m = (ay)£l- =1 is said to be 

(i) stochastic if its elements are non-negative and each row sum is equal to one; 

(ii) substochastic if its elements are non-negative and each row sum is less or equal to one; 

(iii) doubly stochastic if its elements are non-negative and each row and column sums are 



(N 



^1- 



equal to one. 



In [1] G.D Birkhoff characterized the set of extreme doubly stochastic matrices. Namely his 
result states as follows: the set of extreme points of the set ofmxm doubly stochastic matrices 
coincides with the set of all permutations matrices. 

One can consider a generalization of Birkhoff's result in two directions. In the first direction, 
one may consider the description of all extreme points of the set of infinite doubly stochastic 
matrices, and in the second one, one may consider the description of all extreme points of the 
set of nonlinear doubly stochastic operators. 

Concerning the first case, in [10], [13], the Birkhoff's problem have been solved, i.e. it was 
proved that there are no extreme points of the set of all infinite doubly stochastic matrices 
except the permutation matrices. In [18, 19] Yu. Savarov has shown that, under certain 
conditions, Birkhoffs result on doubly stochastic matrices remains valid for countable families 
of discrete probability spaces which have nonempty intersections. Let us also mention some 
other related results. For example, in [14] it was proved that an extreme doubly substochastic 
matrix is a subpermutation matrix. For its generalization to arbitrary marginal vectors see 
[2], for the finite dimensional case and [6, 15], for the infinite dimensional case. In [7], [8] the 
extreme symmetric stochastic and substochastic matrices, respectively, were determined. These 
results were generalized to finite symmetric matrices with given row sums by R. A. Brualdi [2] . 
Finally in [9], [5] the extreme points of the set of infinite symmetric stochastic matrices with 
given row sums were described. 

The present paper is related to the Birkhoff's problem for nonlinear doubly stochastic opera- 
tors. In the this case, we will face with a few contretemps. In fact, first of all, we should define a 
conception of stochasticity for nonlinear operators. We then should define doubly stochasticity 
of nonlinear operators. After all of these, we can consider Birkhoff's problem for nonlinear 
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operators. However, a conception of doubly stochasticity for nonlinear operators can be given 
by different ways. Here, we shall present one of conceptions of doubly stochasticity in nonlinear 
settings introduced in [3]. 

Let us recall some necessary notions and notations. 

Let I m — {1, 2, • • • , m) be a finite set and S m ~ x be an m — 1 dimensional simplex, i.e., 

x = ( Xl ,x 2 ,--- ,x m ) €K m : J2 Xi = 1 > X *-° 

Every element of the simplex S' m_1 can be considered as a probability distribution of the finite 
set I m . Hence, the simplex S*" 1-1 is a set of all probability distributions of the finite set I m . 

Any operator V which maps the simplex S m ~ l into itself is called a stochastic operator. 

For a vector x — (x±, X2, ■ ■ ■ , x m ) G M. m we denote by xj, = (x[q, xp], • • • , X[ m ]) the vector with 
same coordinates, but sorted in non-increasing order X[i] > xp] > ■ • ■ > X[ m ]. For x, y G M. m , we 
say that y is majorized by x (or x majorizes y), and write y -< x if 

k k mm 

^2 y\<\ - S for a11 k = 1? m ~ 1? and = 

i=i i=i «=i i=i 

The Hardy-Littlewood-Polya theorem (see [12]) says that y is majorized by x, i.e., y -< x if and 
only if there exists a doubly stochastic matrix A m such that y = A m x. As a corollary we can 
get that a matrix A m is a doubly stochastic if and only if A m x -< a; for any x G M m . Thus, we 
can give another equivalent definition of the doubly stochasticity of the matrix as follows: a 
matrix A m is called doubly stochastic if A m x -< x for any x G M m . 

Based on this result, in [3] it has been introduced a definition of doubly stochasticity for 
nonlinear operators. Namely, a stochastic operator V : S" Tl ~ 1 — > S 171 ' 1 is called G— doubly 
stochastic if Vx -< x for any x G S m ~ l . An advantage of this definition is that G— doubly 
stochastic operators are well defined for any kind of nonlinear stochastic operators, even though 
the forms of nonlinear operators are not polynomial There is another way to define the notion 
of doubly stochasticity for quadratic stochastic operators. 

Among nonlinear operators, the simplest one is a quadratic one. Such a quadratic operator 
V : lR m — > lR m can be given as follows 

(mm m 
^ ^ A-ij^XiXji ^ ^ Aij^XiXj , • • • , ^ ^ Aij m XiXj 
i,j=l i,j=l i,j=l 

where Ay = (A?,fc)i},fc=i ^ s a cubic matrix. One can see that every quadratic operator is 

uniquely defined by a cubic matrix Ay. In fact, if we denote A™ = (Aj,fc)^=i then the 
quadratic operator has the following form 

Vx={(AMx,x),..- ,(A^x,x)), 

where, (•, •) is the standard inner product in IR n . 

In what follows, we shall use the notation (^Am \ • ■ ■ | A^^j for the quadratic operator V. 

In this paper we attempt to deal with Birkhoff 's problem for quadratic G— doubly stochastic 
operators \ 



1 Here for the sake of completeness we should mention that there is also another way to define quadratic doubly 
stochastic operators is the following sense: a quadratic operator V is called Z— doubly stochastic if its cubic 



matrix Ay satisfies the following conditions ^ij,k — ^-%j,k — ^ij,k — b Mj,k > 0, for all k = l,m. 

i=i ' j=l ' k=l 

One can easily check that if V is a quadratic Z— doubly stochastic then V is stochastic. Note that Z— doubly 
stochasticity of quadratic operators differs from G— doubly stochasticity. However, the disadvantage of the 
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Let us define the following sets 



U m — \ A m — ( y dij) i j = i : ciij — dj 



i > 0, ^2 dij < \a\, Va C I m > , 

i,jea ) 



U m =<jA m GU m : J2 



ciij = m 



where \a\ stands for a number of elements of a set a. 

In [4], the investigation of extreme quadratic G— doubly stochastic operators has been started. 
One of the main results of the paper [4] is that if a quadratic stochastic operator V = 

(^Am | • • • | A^^j is G— doubly stochastic, then the corresponding m x m matrices A$ be- 
long to the set U m for any k = l,m. In other words, the set of all quadratic G— doubly 
stochastic operators is a convex subset of the set XJ m x x U"\ It is clear that the sets U m 

m 

and U" 1 x x are convex. Before studying extreme points of the set of all quadratic 

TO 

G— doubly stochastic operators, it is of independent interest to study geometrical structures 
of U m . The relationship between extreme points of the set of all quadratic G— doubly sto- 
chastic operators and extreme points of U m is given in [4]: let V = [A^ | ••• | A^) be 

a quadratic G— doubly stochastic operator. If any m — 1 matrices of the matrices {Am^}™ =1 
are extreme in the set U m then the corresponding quadratic G— doubly stochastic operator 

V = (^Am | • • • | A^f^j is extreme in the set of all quadratic G— doubly stochastic operators. 

This result encourages us to study extreme points of U m . 

One of the crucial point in the Birkhoff 's problem for quadratic C7— doubly stochastic opera- 
tors is a notion of G— decomposition of symmetric matrices. Namely, let Ai m xm be a set of all 
m x m matrices, and Q C -M mX m be a convex bounded polyhedron. 

Definition 1.1. We say that a matrix A m has a G— decomposition in a class Q if there exists 
a matrix X m G Q such that 

Am = Xm+ 2 Kl , (1) 

By Q s we denote the class of all such kind of matrices A m . The set Q s is called the symmetriza- 
tion ofQ. 

Note that such a notion of G— decomposition is related to certain problems in the control 
theory. 2 One of the fascinating result of the paper [4] is the following one. 



this definition is that Z— doubly stochastic operators are only well defined for polynomial nonlinear stochastic 
operators. In [11] it was concerned with possible generalizations of Birkhoff's problem to higher dimensional 
stochastic matrices and provided lots of criteria for extremity of such matrices. However, the provided criteria 
given in [11] is difficult to check in practice. Therefore, up to now, there is no a full explicitly description of 
extreme higher dimensional stochastic matrices. Particulary, there is not a full explicit description of extreme 
quadratic Z— doubly stochastic operators as well. In [17], it was checked one class of quadratic Z— doubly 
stochastic operators to be extreme. 

2 Recall that the Lyapunov equation has a form Y m X m + X^F m = A m which appears in many branches of 
the control theory, such as stability analysis and optimal control. Here A m ,X m are m x m given matrices and 
Y m is an m x m unknown matrix. If one considers the Lyapunov equation with respect to X m and Y m = \ l m , 
where 1 TO is the unit matrix, then we get (1). Some observations show that if a symmetric matrix A m has the 
decomposition (1) for the special matrix X m then certain problems of convex analysis can be easily solved. 
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Theorem 1.2. [4] Let A m be a symmetric matrix. Then the following statements are equivalent: 

(i) The matrix A rn belongs to the set XJ m ; 

(ii) The matrix A m has a G— decomposition in class of stochastic matrices; 

(iii) The inequality X[ m ] < (A m x,x) < x^ holds for all x G S 171-1 . 

To be fair, we would say that in the paper [4] the provided proof of the part (i) <^ (ii) of 
Theorem 1.2 had some gaps. To clarify and fill those gaps, we aim to write the present paper 
as a complementary one to [4]. Here, we are going to give a complete proof of Theorem 1.2. 
Moreover, we shall generalize it for substochastic matrices as well. As we already mentioned 
above there is a relationship between extreme points of the set U m and the set of quadratic 
G— doubly stochastic operators. However, the extreme points of the set U m were not described 
in [4]. Therefore, we are going to deeply study geometrical and algebraical structures of the 
sets U m and U m . This paper contains many results which are of independent interest. 

Let us briefly explain the organization of the paper. The main results of this paper is the 
following theorem. 

Theorem 1.3. The following statements hold true: 

(i) A symmetric matrix A m has a G— decomposition in the class of stochastic matrices if 
and only if A m belongs to the set XJ m ; 

(ii) A symmetric matrix A m has a G— decomposition in the class of substochastic matrices 
if and only if A m belongs to the set U m . 

The strategy of the proof of the main results is the following: in both cases it is enough to 
prove the assertions of the theorem for extreme points of the sets U m and U m . Then we shall 
employ the Krein-Milman's theorem to prove the theorem in general setting. First of all, we 
shall prove the case (i) for extreme points of U m , then using canonical forms of extreme points 
of U m we reduce the case (ii) to (i). Therefore, our first task is to study extreme points of the 
sets U m and U m . 

In section 2 we shall provide criteria for extreme points of the sets U m and U m . We stress 
here, that the implication (i) =>- (iii) of Theorem 2.16 was stated in the paper [4] without 
any justification. Actually, this implication was a main point of the implication (i) <^ (ii) of 
Theorem 1.2. In this section, we shall justify it, moreover we show that the inverse implication 
(iii) =>- (i) of Theorem 2.16 is also valid. The main results of this section are Theorems 2.14 
and 2.16. 

In section 3 we shall study explicit forms of the extreme points of the sets U m and U m , 
respectively. The results of this section would be used to solve Birkhoff 's problem for quadratic 
G— doubly stochastic operators. The main results of this section are Corollaries 3.11 and 3.12. 

In section 4 we shall study canonical forms of the extreme points of U m which is an extremely 
important to prove the case (ii) of Theorem 1.3. Using the canonical forms of the extreme points 
of U m we are able to reduce the case (ii) of Theorem 1.3 to the case (i) of the same theorem. 
The main results of this section is Corollary 4.9 

In section 5 we shall prove the main results of this paper. They are provided by Theorems 
5.2 and 5.3. There, by means of the results of section 2 we first prove Theorem 5.2 then again 
using the results of section 4 we reduce the proof of Theorem 5.3 to Theorem 5.2. 

2. Some criteria for extreme points of the sets U m and U m 

In this section we want to give some criteria for extreme points of the sets U m and U m . 
Moreover, we provide a proof of some facts which were not proven in [4]. 
It is clear that XJ m is a convex set and U m is a convex subset of U m . 
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One can easily see that 



U x ^ U 2 ^ U 3 ^ • • • ^ U m . 

Here, the inclusion U m should be understood in the way that the matrix with smaller 

order can be extended to larger by letting new entries to be zero. More precisely, the inclusion 
Ufe U m means that if Ak G Ufc then there exists A m G U m such that 



A m — 



Qm—kxk 



Okxm—k 
Qm—kxm~k 



where Q mxn means a m x n matrix with zero entries. 
In the same way, we can get that 

ExtrUi <^ ExtrU 2 ^ ExtrU 3 ^ 



■ ^ ExtrlL 



Particularly, we describe 



here, ExtrU^ denotes the set of the extreme points of U^. 

We are going to study a geometrical structure of the set U m 
extreme points of U m . Let us recall some well-known notations. 

A submatrix A a of A m = (ajj)™- =1 is said to be a principal submatrix if A a 
all entries indexes of A a belong to a(c I m ). 

The proof of the following proposition is straightforward. 

Proposition 2.1. The following statements hold true: 

(i) If A m G U m , then its any principal submatrix of order k belong to Ufc, 1 < k < m; 



i.e. 



ii 



If a principal submatrix A a of A m is extreme in U\ a \, for some a C I m , then for any 
matrices A' m , A" m G U m satisfying 2A m = A' m + A" m , one has A a = A' a — A" a . 



Proposition 2.2. A matrix A rn belongs to ExtrU m if and only if for any entry of A m , 
there exists a G I m such that is an entry of the principal submatrix A a with A a G ExtrU| Q |. 

Proof. If part. Let us assume that A m is not extreme, that is 2A m = A' rn + A'^, for some 
A' m , A" m G U m with A m 7^ A' m , A m ^ A" m . From the former, we conclude that there is an 
entry aj 0j0 such that di j ^ a i j - According to the condition, there exists an extremal principal 



submatrix A ao in U| Q0 |, containing a,i j . Then due to Proposition 2.1 we get A ao 



A' 



a! in which is a contradiction. 



and hence a io j 

Only If part. Suppose that A. 
A a G ExtrU| Q |, for any entry 



A" 



G ExtrU m . By putting a = I m we get a 



A a and 
□ 



Remark 2.3. Note that the provided criterion is somehow difficult to apply in practice. The 
reason is that sometimes a given matrix may not have extreme proper principal submatrices. 
Let us consider the following m x m matrix 



i 

2 


1 

2 




1 

2 








1 

2 



V 




1 

2 






























\ 



1 

2 







1 

2 



1 

2 





/ 



For this matrix, the problem of finding of its extreme proper principal submatrix coincides with 
the problem of showing its extremity. Further, one can show that the matrix N m is not extreme. 

However, there are some benefits of the provided criterion in terms of studying some prop- 
erties of extreme matrices. The following corollaries directly follow from Proposition 2.2. 
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Corollary 2.4. A matrix A m is not extreme in XJ m if and only if there exists an entry a io j 
such that any principal submatrix A a containing this entry, is not extreme in U\ a \. 

Corollary 2.5. If every 2x2 principal submatrix of a matrix A m is extreme in U2 then the 
matrix A m itself is extreme in U m . 

Remark 2.6. The converse of Corollary 2.5 is not true. For instance, the matrix 

M 3 

V i 1 J 

! 
1 

2 

Let us now present some facts. We are not going to prove them because of their evidence. 
Proposition 2.7. The following assertions hold true: 




is extreme in U 3 , however it has a 2 x 2 non extreme principal submatrix ( 1 ^ 



(i) If A m G U m , then < < 1 for all i,j = l,m; 

(ii) Let matrices A m , A' m , A" m G U m satisfy 2A m = A' m + A" m . If a iojo = 1 V for some 
io,jo £ Im-, then ai j = a- M0 = a" M0 - Here and henceforth a — feV c means either a = b 
or a = c; 

(hi) Any matrix A rn G U m having entries being equal to either 1 or is extreme in U m . 

Let us introduce the following useful conception. 

Definition 2.8. Let A m G U m . An index set a C J m sazd to &e saturated, for the matrix 
A rn , whenever Yl a ij = \ a \- A principal submatrix A a , corresponding to the saturated index 

i,j€a 

set a, is called a saturated principal submatrix. 

An advantage of the given conception is that using induction with respect to the number of 
saturated index sets we can easily prove lots of properties of U m . Moreover, it is an appropriate 
conception to formulate some facts regarding extreme matrices of U m . 

Proposition 2.9. Let a, j3 C I m be saturated index sets of A rn G U m . Then the following 
assertions hold true: 

(i) If a fl j3 7^ 0, then ad f3 is a saturated index set for A m ; 

(ii) a U j3 is a saturated index set for A m . 

Proof. Let 7 = a fl f3, for any index sets a, (3 C I m . Then, one can easily check the following 
equality 

a ij = a ij + a «i + a ij + a u — a «r (2) 

ij'GaU/3 iec«\7 »e/3\7 *j£/3 *j€7 

j'6/3\7 j6a\7 

Now suppose that a, (3 C I m are saturated index sets of a matrix A m G U m and 7 7^ 0. We 
only consider a case when 7 7^ a, 7 7^ /3 otherwise the theorem is evident. It is clear that 
^2 a ij — Hence by means of (2) we have 

\a U p\ = \a\ + \(3\ - |7| < ai i + Y a ^ ~ a ^ 

ij'e« i,jep ij'67 

iJeaU/3 
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This yields that 

i,j£aU/3 

Therefore a fl j3 and a U j3 are saturated index sets for A m e U m . □ 

According to Proposition 2.9, a class of all the saturated index sets of a given matrix is 
closed with respect to the operations of union and intersection. That is why, this property of 
the saturated index sets implies a reason to introduce the following 

Definition 2.10. A saturated principal submatrix of a matrix A m e U m containing an entry 
aij is called saturated neighborhood o/ay, and the order of such a saturated principal submatrix 
is said to be its radius. Saturated neighborhoods of aij with minimal and maximal radiuses are 
called a minimal and a maximal saturated neighborhoods of a^, respectively. 

Remark 2.11. If an entry of A rn e U m has a minimal or a maximal saturated neigh- 
borhoods, then they are uniquely defined. Let us show uniqueness of the minimal saturated 
neighborhood of a^. Assume that there are two minimal saturated neighborhoods A a , A a i of a^, 
and we denote the corresponding saturated index sets by a and a' . Since aC\ a' ^ and A a , 
A a i are minimal saturated neighborhoods, due to Proposition 2.9 (i), a Ha' is a saturated index 
set, and \a\ < \a fl a'\, \a'\ < \a fl a'\. Therefore, a — aHa' — a'. 

Using the same argument with Proposition 2.9 (ii), one can get the uniqueness of the maximal 
saturated neighborhood of aij. 

We would like to emphasize that the minimal saturated neighborhood plays an important 
role, for geometrical structures of the set U m whereas the maximal saturated neighborhood 
plays as crucial point for its algebraical structures. 

If an entry a^ has a saturated neighborhood then, since the matrix A m e U m is symmetric, 
an entry has also the same saturated neighborhood. That is why, henceforth, we only 
consider saturated neighborhoods of entries a^ in which % < j. 

Let us observe the following: assume that A a is a principal submatrix of a matrix A m e U m 
and entry of A a . Then, in general, the minimal saturated neighborhood of a^ in 

the matrix A rn does not coincide with its minimal saturated neighborhood in the principal 
submatrix A a . For this, one of the main reasons is that the entry a^ may have a minimal 
saturated neighborhood in A m , but may not so in A a . We can see this picture in the following 
example: let 

/o|oooo\ 

\ | 



A,= 



| | 

\ 1 



1 

\ 

1 n 1 



If we consider the principal submatrix A ao — ( 2 ^ 2 I ' wnere a o — {1)2,3}, then the 

§ 



element 023 = \ does not have a minimal saturated neighborhood in A ao , but it has the minimal 

\ I in the given matrix A 6 

6 h 



i 

saturated neighborhood Ap — I \ \ I in the given matrix A 6 , where (3 = {2,3,4}. 
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Fortunately, if a principal submatrix A a of a matrix A m is saturated, then a minimal saturated 
neighborhood of any element a^, where i,j G a, in the matrix A m , coincides with its minimal 
saturated neighborhood in the principal submatrix A a . Namely, we have the following 

Lemma 2.12. Let A a be a saturated principal submatrix of a matrix A m G U m and be 
any entry of A a . Then there exists a minimal saturated neighborhood of in A m , and it is a 
principal submatrix of A a . Moreover, it is a minimal saturated neighborhood of in A a . 

Proof. Suppose that A a is a saturated principal submatrix of A m G U m and is an entry of 
A a . Then has at least one saturated neighborhood in A m , which is A a . Therefore, it has a 
minimal saturated neighborhood in A m . By a' we denote a saturated index set corresponding 
to its minimal saturated neighborhood in A m . We must show that a' C a. Indeed, since 
aflft'D {i, j} and both a and a' are saturated index sets, then according to Proposition 2.9 
the set a D a' is saturated as well. Further, since a' corresponds to the minimal saturated 
neighborhood of in A m , it follows that a'Ca'HaCaas desired. Moreover, it immediately 
follows from the definition of the minimal saturated neighborhood and a' C a that the minimal 
saturated neighborhood of a^-, in A m , is its minimal saturated neighborhood in A a . □ 

It is worth mentioning that due to Proposition 2.7 we shall deal with entries such that 

< CLij < 1. 

Let us present some criteria for the extremity of a matrix A m G U m . 

Theorem 2.13. Let A m be an element of\J m . Then the following assertions are equivalent: 

(i) The matrix A rn is an extreme point of\J m ; 

(ii) Every entry of A rn with < < 1 has at least one saturated neighborhood, and 
minimal saturated neighborhoods of any two entries a^, a^y with < 0^,0^' < 1 do 
not coincide. 

Proof. (i)=> (ii). Let A m e ExtrU m . Let us prove that every entry of A m with < a^- < 1 
has at least one saturated neighborhood. We suppose the contrary, i.e., there exist < a io j < 1 
having no saturated neighborhoods, which means, for any a G I m with a D {io,jo} one has 



We know that there are two cases either i 7^ j or i = j . Let us consider the case i 7^ j , 
the second case can be proceeded by the same argument. 

Since I m is a finite set, then a number of inequalities in (3) is finite. That is why, there exists 
< e < 1 such that 




(3) 



i,j€a 



< a iojo + e < 1 



< Oj j Q — Eq < 1 , 





i,jea 



for all a D {i , j }. 
Define two matrices 




and A" 




as follows 
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Due to the choice of e , we get A' m , A" m e U m and 2A m = A' m + A" m which refutes the extremity 
of A m . 

Let us show that, for any two entries a^, a^y with < ai^ayy < 1, their corresponding 
minimal saturated neighborhoods do not coincide. We suppose the contrary, i.e., for two entries 
< a iQ j Q < 1 and < o,i' j> < 1 their corresponding minimal saturated neighborhoods coincide. 
This means that there exists «o C I m with {i , j , i' ,f } C «o such that 



a 



|«o| 



We suppose that i ^ j , i' Q ^ j' . For the other possible cases one can use the similar 
argument. 

Since A ao is a common minimal saturated neighborhood of entries a io j and at y , then due 
to the minimality of «o it is clear that, for any a with either a D {io,jo} or a D {i' , f } and 



'J 



we have « C a. In other words, any saturated neighborhoods of a io j or contain both of 
them. 

Let us consider all index sets j3 D {^o, jo} and D {i )Jo} sucn that 



< 



< 



(4) 



Since I m is a finite set then a number of inequalities in (4) is finite. Therefore, one can find 
< Eq < 1 such that 



< at 



"iojo ± £ o < 1, < ±e < 1, 

i,j6/8 ij'6/3' 

for all (3 D {i , jo} and /?' D {i , f } satisfying inequalities (4). 

Define two matrices A' m = (a^)™ ?=1 and A" rn = (a^)™ =1 , as follows 



""WO 



+ ^0, 

°io*o ~~ £ 0) 
a i' f + £ 0, 



{hJ) = Vo,3o) 
(ij) = (jo,io) 
= (i'oJ'o) 

= (*o, Jo) 

= (io,«o) 

= (*o,Jo) 

^ (io,jo) V (j ,i ) V (z , j ) V (j ,i ) 



Let us show that A' m , A" m e U m . In fact, the matrices A' m , A'^ are symmetric. We shall show 
that A' m G U m . By the same argument, one can show that A" m e U m . 

Let /3 be any subset of I m . We want to estimate the sum a!^. We are going to consider 

the following cases. 

Case I. Let (5 D {i , j } and f3 D {i' , j' }. Then one gets a'^ = < 
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Case II. Let f3 £ {i ,j } and f3 £ {i' ,f }. Then £ <■ = E «u < 

Case III. Let ^ {i„, jo} and (3 D {i' ,j' }. Then' £ e^- =' £ a^- - 2e < 

Case IV. Let /3 D {i , jo} and /3 ^ {i >.7o}- Tnen * E 4? = E + 2e o- It is clear that 

i,j£/3 i,j<E0 

the set /3 could not be a saturated index in the matrix A m . In fact, every saturated index in 
the matrix A m containing the set {io,jo} should contain the set {ib,f }- However, j3 {i'o,j'o}- 
Therefore, due to choice of e we have E a'^ = E a ij + ^ £ o < 

Consequently, the constructed matrices A' m , A" m belong to U m , and 2A m = A' m + A'^ which 
refutes the extremity of A m . 

(ii)=>- (i). Suppose that the assertions (ii) are satisfied. Let us prove A, m e ExtrU m . We are 
going to show it by using induction with respect to the order of the matrix A m . 

Let m = 2. If the entries of the matrix A2 G U2 are either 1 or 0, then according to Proposition 
2.7 (iii) A 2 is extreme. Suppose that there exist at least one entry of A 2 with < < 1. 
Then it is obvious that such kind of entries' saturated neighborhood's radius is greater or equal 
to 2. Consequently, A 2 is a minimal saturated neighborhood for all < < 1. It follows from 
(ii) that there is only one entry < < 1 (of course, we only consider entries with i < j) 
and the rest entries are either 1 or 0. After small algebraic manipulations, we make sure that 
A 2 is an extreme matrix in U2. 

We suppose that the assumption of the theorem is true for all m < k — 1, and we prove it 
for m = k. 

If the entries of the matrix A k e U k are either 1 or then according to Proposition 2.7 (iii) 
A k is extreme. Suppose that there exist some entries < < 1. It follows from (i) that every 
entry < < 1 has a minimal saturated neighborhood, and we denote it by A a ^ a ..^, its radius 
by r{a ij ). 

Let us consider entries with < < 1, r(a^) < k—1. Since A a ( a ..) is a saturated principal 
submatrix of A m , then according to Lemma 2.12, the minimal saturated neighborhoods -A Q ( 0j , ,) 
of entries < ayy < 1 in A a ^ a ..^ are principal submatrices of A a ( a ..y So, theorem assertions 
(i),(ii) are satisfied, for A a ( aij ), and order of such a matrix A a ( a .^ is less or equal to k — 1. 
Hence, according to the assumption of the induction we obtain that A*(a i3 ) e ExtrU| Q ( aij )|. So, 
minimal saturated neighborhoods A a ^ a ..) of entries < < 1 of A k with r(ay) < k — 1 are 
extreme in U| a ( aiJ .)|. 

Let us consider entries with < < 1, r(ajj) = k. It follows from (ii) that there is only 
one such kind of entries. 

Suppose that the matrix can be decomposed as 2A k = A' k + A k , where A' k ,A k e U k . If the 
entries of A k are either 1 or then according to Proposition 2.7 (ii) we get 
Further, for the entries < < 1 with r(dij) < k — 1, since their minimal saturated neighbor- 
hoods A a ( ai] ) are extreme in U| a ( 0jJ .)|, according to Lemma 2.1 we have A a{aij ) = A' a ^ = A'^ a ^, 
particularly, = a'^ = a'-j. Now, we must to show = a'^ = a'-j for an entry < < 1 
with r(ciij) = k. We already mentioned that there is only one such kind of entries, we denote 
it by a io j and its minimal saturated neighborhood is A k . Since 2A k = A' k + A' k then A' k and 
A' k are also saturated matrices. We already know that for entries of A k with 7^ (io,jo) 
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and ^ (j ,i ), one has Therefore, we get 

k k 

= k-k=^Ta ij -Y^ a 'ij = 2 ( a wo - a< j ), 

k k 

= k — k = ^ ciij — ^ a"j = 2(a io j — a" ojo ). 

Consequently, a io j = a- oJo = a" j - All these facts bring to a conclusion that A k is a extreme 
matrix in U^. □ 

Theorem 2.14. Let A m be an element of\J m . Then the following conditions are equivalent 

(i) The matrix A m is an extreme point o/U m ; 

(ii) Every entry of A m with < < 1 has at least one saturated neighborhood and its 
minimal saturated neighborhood A a is extreme in Ui a i; 

(iii) Every entry of A rn with < < 1 has at least one saturated neighborhood and any 
its saturated principal submatrix A a is extreme in U\ a \. 

Proof. The implication (iii)=>- (ii) is obvious. Consider the implication (ii)=> (i). Assume 
that A m has the following decomposition 2A rn = A' m + A" m with A' m , A" m e U m . We know 
that if ay = 1 V then Proposition 2.7 (ii) yields If < ciij < 1, since its 

minimal saturated neighborhood A a is extreme in U| a |, then due to Proposition 2.1 one gets 
A a = A' a = A'a, particularly, = a'^ = a"j. These mean A m = A' m = A n m , i.e. A m is a extreme 
matrix in U m . 

(i)=>- (ii). Let A rn e ExtrU m . Then according to Theorem 2.13 every entry with < < 
1 of A m has at least one saturated neighborhood. We want to show that a minimal saturated 
neighborhood A a of such an entry is extreme in Ui q i . Since A a is a saturated principal submatrix 
of A m then using Lemma 2.12 we deduce that, the minimal saturated neighborhood in A m of 
an entry a^y of A a with < a^y < 1 coincides with its minimal saturated neighborhood in 
A a . Therefore, by applying Theorem 2.13 for the matrix A a we conclude that A a e ExtrU Q . 

(i)=> (hi). Let A m e ExtrU m . Then according to Theorem 2.13 every entry a^- of A m 
with < ciij < 1 has at least one saturated neighborhood. We want to show that any its 
saturated principal submatrix A a is extreme in U| a |. Let Ap be a saturated principal submatrix 
of A m G U m . Since Ap is a saturated principal submatrix then using Lemma 2.12 we can 
conclude that, for every entry < ayy < 1 of Ap, its minimal saturated neighborhood in A rn 
coincides with its minimal saturated neighborhood in Ap. Therefore, if we apply (ii)=>(i) to the 
submatrix Ap we get Ap e ExtrU/3. □ 

We are going to describe all extreme points of U m . It is clear that the set U m is a set of all 
the saturated matrices of the set U m . 

Proposition 2.15. A matrix A m e U m is extreme in XJ m if and only if A m is extreme in U m . 
Namely, one has 

ExtrU m = U m n ExtrU m . 

The proof of this proposition is straightforward. 

From Theorem 2.14 and Proposition 2.15 we immediately get 

Theorem 2.16. Let A m be an element of\J m . Then the following conditions are equivalent 

(i) The matrix A m is an extreme point of\J m ; 

(ii) Every minimal saturated neighborhood A a of A m is extreme in U' a '; 

(iii) Every saturated principal submatrix A a of A m is extreme in U' a L 
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It is worth mentioning that in the paper [4] the part (i) =>■ (Hi) of Theorem 2.16 was stated 
without any justification. 

3. Some properties of extreme points of the sets U m and U m 

In this section, we are going to study some properties of extreme matrices of the sets U m 
and U m . The results of this section will be used to solve Birkhoff's problem for quadratic 
G— doubly stochastic operators. 

Let us introduce the following sets 

:= {A m 6U ra | a !J= 0Vl Vi,j G J m }, 
Um 2 ' := ^A m G U m | an = 0, = V ^ Vi, j G I n 

U^' 1 * 1 ' := !^A m eU m \ a« = 0Vl, a - = V i V 1 Vi,jel„ 

Remark 3.1. T/ie following assertions are evident: 

(i) Ufcu? and U^CUS^; 

(ii) If A £ Um 2 ' 1 ' £nen a^- is an integer, for any a C J m . 

Theorem 3.2. 7/ A m G ExtrU m £nen = V 1 and a^- = V | V 1 /or any « ^ j and i, j <E I m . 
In other words, we have ExtrU m C Um 2 ' 1 '. 

Proof. It is enough to show that an and 2ajj are integers, for any i ^ j and i,j G I m . We prove 
it by using induction with respect to the order of A m . 

Let m — 2. Since A2 G ExtrU2, then according to Theorem 2.13, there is at most one entry 
aij with < aij < 1. If there is no such an entry then the claim is obvious. Assume that 
there exists an entry a^ with < a^ < 1. Then the matrix A 2 should be saturated and from 
an + 2ai2 + 022 = 2 we deduce that an, 2ai2, a 2 2 are integers. 

Now suppose that the assertion of the theorem is true, for all matrices A rn G ExtrU„ t of 
order m < k — 1. We prove it for matrices A m G ExtrU m of order m — k. 

Since Aj, G ExtrU m then, due to Theorem 2.14, every entry < a^ < 1 of Aj, has a 
minimal saturated neighborhood A a which is extreme in U| Q |. So, according to the assumption 
of induction, for those minimal saturated neighborhoods with radius less or equal to k — 1 their 
entries an and 2ajj are integers, for any % 7^ j and i,j G Ik- In other words, the assumption of 
induction allow us to say that all an and 2ay (i < j) of entries of A^ are integers except which 
has a minimal saturated neighborhood with radius k. 

Now, assume that there is an entry a^ with < a^ < 1 which has saturated neighborhoods 
with radius equal to k. Then according to Theorem 2.13, there is only one such an entry, we 
denote it by ai j , and the matrix A^ should be saturated. In this case, we already know that 
for entries a^ of A k with 7^ (io,jo), if i — j then an is an integer, and if i 7^ j then 2a^ is 
an integer. Then from 

k k 

} ] a ij = ^ ] a ii + 2 ^ ^ Ojj = k 

i,j=l i=l i<j 

we conclude that if io = jo then ai i is an integer, and if io 7^ jo then 2a,i j is an integer as 
well. All of these mean that an and 2a^ are integers, for any i 7^ j and i,j G Ik- □ 
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Now, combining Proposition 2.7 (iii) with Theorem 3.2 we have 

(oil) 



U^' 1 ) C ExtrU m C IC 5 
Remark 3.3. If m = 2 then it is easy to show the following equalities: 

(i) n ExtrU 2 = 02x2; 

(ii) U^ |) UExtrU 2 = U^^' 1) ; 

here, as before, 02x2 ^2x2 zero matrix. 



(5) 



Our next aim is that Remark 3.3 (i) holds true for any m. To this end, we introduce the 
following useful conception. 



Definition 3.4. A matrix A m G U 



(o,ii) 



is said to be F m — matrix if A m is saturated, and having 



no other saturated principal submatrices. 



(o,|) 



Proposition 3.5. Let A m be an F m — matrix. If m > 3 then A m G U 

Proof. Suppose that A m is an F m — matrix and m > 3. If a — {i} then it follows from the 
definition of the F m — matrix and Remark 3.1 (ii) that an < 0, which means an = for all 
i G I m . If a = {i,j} with i ^ j, then with the same reason as the previous case, from the 
inequality 

^ ^ Q>ij = Q>u ~\~ Oijj -\- 2aij ^ 1, 



we find out that a^ — V for alH ^ j and i, j G I m . This means that A m G U 
Remark 3.6. For m = 2, there is only one F 2 — matrix which is 



(0.3) 



□ 



From now on, we shall only consider the case m > 3. 
Example 3.7. Let us consider the following m x m matrix 

( \ ••• 





i 

2 





1 

2 


1 

2 





1 

2 










V 




1 

2 




























1 

2 





1 

2 







1 

2 



1 

2 





It is easy to see that if m > 3 then N m is an F m — matrix. 



Proposition 3.8. Let A m be an F m — matrix. If m > 3 then A m is not an extreme element of 
U m . 

Proof. Suppose that A m is an F m — matrix and m > 3. Then, due to Proposition 3.5, we have 
(o -) 

A rn G U m ' 2 . By the definition of the F m — matrix, the number of entries of A m , which are 1/2, 
is equal to m (of course, we are only speaking about such entries in which i < j). On the 
other hand, a minimal saturated neighborhood of any nonzero entry is A m itself. So, there are 
at least two nonzero entries whose minimal saturated neighborhoods coincide. Then Theorem 
2.13 brings to a conclusion that the matrix A m is not extreme in U m . □ 

r(o,i,i) 



Proposition 3.9. Let A m G U r 

A m is not extreme in U m . 



If A m has a principal F k — submatrix, where k > 3, then 
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Proof. Let A m G uJ!' 2 ' 1 ' and be a principal Fk~ submatrix of A m . Since k > 3 then, due 

(0 A) 

to Propositions 3.5 and 3.8, we get A k G U fc " 2 and A fe ^ ExtrU^. Then a number of entries 

of Ak, which are equal to -, is equal to k and a minimal saturated neighborhood of such an 

entry is Ak, which is not extreme in Then, according to Theorem 2.14, we conclude that 
the matrix A m is not extreme in U m . □ 

Now, we are ready to formulate one of important properties of the extremal matrices in U m . 
Theorem 3.10. The following equality is satisfied for any m 

ExtrU m n U^' 2 "' = Q mxm , (6) 

here, as before, Q m xm is zero matrix. 

Proof. Let m > 3, otherwise Remark 3.3 (i) yields the assertion. Suppose that A m G ExtrU m 

(0,1) (o,-) 
and A m ^ Q m . We want to show that A rn f. U m ' 2 . We assume that A m G Um' 2 • Since A m 

(o -) 

is a nonzero extreme matrix in U m ' 2 , then due to Theorem 2.14, there exists at least one 
entry = | and its minimal saturated neighborhood A a should be extreme in U Q . On the 

(0,1) 

other hand, since A m G U m ' 2 then the minimal saturated neighborhood A a of that entry 
(i.e. aij = |) is F\ a \— matrix which is not extreme U a . This contradiction proves the desired 
assertion. □ 



Corollary 3.11. Let A m G ExtrU m , then the following assertions hold true: 

(i) One has A m G Um 2 ' 1 '; 

(ii) Every entry = \ has at least one saturated neighborhood; 

(o -) 

(iii) For any saturated principal submatrix A a , one has A a f. U| a j 2 . 

By using Corollary 3.11 and Proposition 2.15, we get the following description of the extreme 
points of U m . 

Corollary 3.12. Let A m G ExtrU™. Then the following assertions hold true: 

(i) One has A m G Um 2 ' 1 '; 

(o -) 

(ii) For any saturated principal submatrix A a , one has A a U| a j 2 . 
It seems the following conjectures hold true. 



Conjecture 3.13. Let A m G ExtrU m if and only if the following assertions hold true: 

(i) One has A m G U^' 1 '; 

(ii) Every entry aij = \ has at least one saturated neighborhood; 

(o -) 

(iii) For any saturated principal submatrix A a one has A a U| a j 2 . 

Conjecture 3.14. Let A rn G U m . Then A m G ExtrU™ if and only if the following assertions 
hold true: 

(i) One has A rn G U^' 1 '; 

(o -) 

(ii) For any saturated principal submatrix A a , one has A a f. U| a j 2 . 
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Proof, (i). Due to Corollary 3.11 (i) and (ii), we have that A m G U m ' 2 ' and every entry 



4. Canonical forms of extreme points of the sets U m and U m 

In this section, we are going to describe location of nonzero entries of extreme matrices and 
canonical forms of extreme matrices of the sets U m and U m . Based on the canonical forms of 
extreme points we are going to study an algebraic structure of the sets U m and U m (see sec. 
5). 

Proposition 4.1. Let A m G ExtrU m . Then the following statements hold true: 

(i) Every nonzero entry of A rn has a minimal saturated neighborhood; 

(ii) Every nonzero entry of A m has a unique common maximal saturated neighborhood. 

r(°>i!) „ _ I 

2 

has a minimal saturated neighborhood. Now, let us consider such entries with = 1. Then, 
there are two cases either % = j or i ^ j. In both submatrix A a with a = {i,j} of A m 

is a minimal saturated neighborhood of a^- = 1. 

(ii). Due to (i) every nonzero entry of A m has a minimal saturated neighborhood. By 
«i, c*2, • • • , otk we denote saturated index sets, corresponding to these minimal saturated neigh- 

k 

borhoods of nonzero entries of A m . According to Proposition 2.9 (ii), an index set a = [J a>i 

i=i 

is saturated. If we consider a principal submatrix A a of A m , corresponding to the saturated 
index set a, then A a is saturated. Since A a contains all nonzero entries of A m , therefore, A a is 
a common maximal saturated neighborhood of every nonzero entry of A rn . The uniqueness is 
trivial. □ 

Corollary 4.2. If A m G ExtrU m; then there exist two index sets a and a' such that aUa' = I m , 
a fl a' = 0, and satisfying the following conditions: 

(i) A a is a saturated principal submatrix of A m ; 

(ii) A a contains all nonzero entries of A m ; 

(iii) A a i = Q\ a '\ and a^i = a^j = for any i,j G a, and i',f G a'. 

Lemma 4.3. Let A m G ExtrU m . A matrix A m is saturated if and only if every row of A m has 
at least one nonzero entry. 

Proof. If part. Let A rn G ExtrU m , and its every row has at least one nonzero entry. Then, due 
to Corollary 4.2 (i) and (ii), there exist two index sets a and a' with aUa' = / m , a fl a' = 
such that A a is a saturated principal submatrix of A m containing all nonzero entries of A m . 
Moreover, since every row of A m has at least one nonzero entry, it follows form Corollary 4.2 
(iii) that a' = 0, which means A m = A a is a saturated matrix. 

Only if part. Let A m be a saturated extreme matrix in U m . We must to show that every row 
of A m has at least one nonzero entry. Assume the contrary i.e. there is a row of A m with zero 
entries. We denote it by iq. In this case, symmetricity of A m implies that ig h column of A m has 
zero entries as well. Let ao = I m \ {io}, then from A m G U m , one gets 



E 



dij < \oio\ — m — 1. 
On the other hand, since A m is saturated, it follows that 

m 

i,j=l ijeao j£ao i,j£ao 

This contradiction shows that every row of A m has at least one nonzero entry. □ 
By using Proposition 2.15 and Lemma 4.3 one can get the following 
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Corollary 4.4. Let A m G ExtrU m . Then the matrix A m is extreme in XJ m if and only if its 
every row has at least one nonzero entry. 

m 

Corollary 4.5. Let A rn be a nonsaturated extreme matrix in XJ m and K := m— Yl a %i- Then, 

there exist two index set a and a' with \a'\ = K, a = I m \a', satisfying the following conditions: 

(i) A a is a saturated principal submatrix of A m ; 

(ii) A a contains all nonzero entries of A m ; 

(iii) Every i row and i column of A m consists zeroes, for any i G a'; 

We will introduce the following notation. Let A m be a matrix and it be a permutation of the 
set I m — {1, 2, • • • , m}. Define a matrix as follows 

^7r(m) = ( a 'ij)™j=n a 'ij = a ir(i)Tr(j) V i,j = 1,771. 

Proposition 4.6. Let A m be a matrix and n be a permutation of the set I m = {1,2, •• • , m}. 
Then the following assertions hold true: 

(i) if A m G U m then A n{m) G 

(ii) if A m G ExtrU m then A^ m) G ExtrU m ; 

(iii) if A m G U m then A n(m) G U m ; 

(iv) if A m G ExtrU" 1 then A n{m) G ExtrU" 1 . 



Proof, (i). Let A m G U m . Due to — aj i: V i,j — l,m we have 

/ I y • • 1 

O jj ^7r(i)7r(jr) ^7r(jr')7r(j) a jri; » 2, J 1,777. 

This means that (A,.( m ))* = A^.^. Let ct C J m , and put /3 = 7r(a). It is clear that = |a|. 
From A rn G U m it follows that 

a>i i = a T(0T(j) = - 1^1 = l a l' 

which implies A n ( m ) G U m . 

(ii) . Let A m G ExtrU m . We want to show that A T ( m ) G ExtrU m . We suppose the contrary, 
i.e., there exist matrices A' m , A" m G U m such that 2A n ( m ) = A' m + A" m and A n ( m ) ^ A' m , 
A n ( m ) 7^ A" m . Let us consider the matrices A' n -i^, A" n -i^ then one has 

2^4 m — ^4 ,r-i( m ) + A 7r-i( m ), j4„.-i( m ) 7^ A m , j4„— i( m ) 7^ A m . 

This contradicts to A rn G ExtrU m . 

m 

(iii) . Let A m G U m , then from = m one finds 

mm, m 
*J = 1 ij'=l i,j=l 

which means A n ( m j G U m . 

Proposition 2.15 with (ii),(iii) yields the assertion (iv). □ 

By means of Proposition 4.6 and Corollary 4.5 we are going to provide a canonical form of 
extreme points of U m . 

m 

Let A m G U m be an extreme matrix in U m and /c = Yl a ij- Then there exists a permutation 
7r of the set I m such that the matrix A n ^ m ) has the following form 

a I Ak ©fcxm-fc 1 lr->\ 

Ar(m) = q q > I 7 ) 

\ ^Jm~kxk ^Jm—kxm-k I 
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here, as before, Qkxm-k, Qm-kxk, and G m _kxm-k are zero matrices and A k is an extreme 
saturated matrix in U fe , i.e., A k G ExtrU fc . The form (7) is called a canonical form of the 
extreme matrix A m . 

In the sequel, without loss of generality, we will assume that an extreme matrix A m has a 
canonical form (7). 

Let a and j3 be two nonempty disjoint partitions of the set I m , i.e., afl/3 = and aU f3 = I m . 
Let A a = (aij)ij ea and Bp = {b%j)i,jep be two matrices. Define a matrix C aU p = (cij) aU p as 
follows: 

dij i,j G a 

• / i a, j 0- 1 1 ' 

% g /?, j G a. 

Proposition 4.7. Let a, /3 fre too nonempty disjoint partitions of I m and A a = (a>ij)ijea, 
Bp = (bij)i,je/3 be two matrices. Let C aU p = (cij) aU p be a matrix defined by (8). Then the 
following assertions hold true: 

(i) if A a G U| a | and Bp G U\p\ then C aU p G U^.^; 
(ii) if A a G ExtrU| a | and Bp G ExtrU^i then C aU p G ExtrU| a |+|0|; 
(Hi) if A a G Ul°l and Bp G TJW then C aU p G TjW+W; 
(iv) if A a G ExtrUl"! and Bp G ExtrU^ then C aU p G ExtrU^+^l. 

Proof. According to Proposition 4.6 we may assume that a — {1, 2, • • • , i} and /3 — {i + 1, % + 
2, • • • , m). Then the matrix Co, U( g = (cij) aU p given by (8) has the following form 

c *» - ( e£ w V ) • (9) 

(i) . Let A a G U| a | and Bp G U^. We want to show that C aU p G UiqI+i^. Indeed, it follows 
from (9) that (C aU pY = C aU p. Let 7 C I m (— a U j3) be any subset of J m . Let 7 Q = 7 fl a, 
7/3 = 7^/3 then 7^ fl ^p = 0, 7 Q U 73 = 7. It is clear that 

This means that C Q u/3 G U|q.|+|j(3|- 

(ii) . Let A a G ExtrU^i and Bp G ExtrU^. We suppose that there exist C' u/3 and C" u/3 
such that 2C aU p = C' aU/3 + C^ u/5 . We may assume that the matrices C' a U(3 and C^ u/3 have the 
following form 

°« °|a|x|/3| I r" - I a HxlP 

U |/3|x|a| U /3 j V U \f>\x\a\ U fi 



then it follows from (9) that 



20, 


a|x|/3| 


~ C \a\x\f3\+ C \a\x\f3 


20, 


/3|x|a| 


= C \f3\x\a\+ C \f3\x\a 




2A a 


= C a + C a , 




2Bp 


— Cp + Cp. 
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Hence, we obtain C' Hm = C^ m = Q\ a \x\p\ and C^ |x|a| = C^ x|q| = Q Wx \ a \. From A a e 
ExtrU w , Bp e ExtrU^ we find C a = C a = A a , Cp = Cp = Bp. This means that C aUj3 — 
C a ui3 = Cau/3, i-e, C Q u/3 £ ExtrU^i-i^i . 

(hi). Let A a e U |a| and Bp e U I/3 L One can see that 



m 



2J=1 i,jea i,je/3 i,jea i,jefi 

this means that C aU p E U |q|+I/3 L 

The assertion (iv) immediately follows from Proposition 2.15 and assertions (ii) and (hi). □ 

Now we are going to consider an extension problem: let A m e U m be a non- saturated matrix. 
Is there a saturated matrix A m+1 e U m+1 containing a matrix A rn as a principal sub-matrix? 
In other words, is it possible to make a non-saturated matrix as a saturated matrix by increasing 
its order? If the extension problem has a positive answer, for the sets U m and U m+1 , then we 
use the following natation U m ^ XJ m+1 . We shall solve this extension problem in a general 
setting. 

Proposition 4.8. Let A m e U m . Then there exists a saturated matrix A rn+ i e U m+1 contain- 
ing a matrix A m as a principal sub-matrix, i.e., U m > XJ m+1 . 

Proof. We shall prove the assertion in two steps. 

Step I. Let us prove that if A m e ExtrU m then there exists a saturated matrix A m+1 e 
ExtrU m+1 containing a matrix A m as a principal sub-matrix, i.e., ExtrU m ^ ExtrU m+1 . 
Indeed, we suppose that the matrix A m e ExtrU m has the following form 

Afc Qkxm—k 
Qm—kxk Sm-hm-f: 

where A fe e ExtrU fe and = k < m. Let us consider the following matrix 



i-m+l-k 



Om—kxm—k (2) 
\2Jlxm-k 



m—kxl 



where (^), , = ( —, ■ ■ ■ ,-] and (^) . . = (^)* , . It is clear that the matrix A m+ i-k 

\2/lxm—k \2 2/ V2/m— fcxl \2/lxm—k m+± re 

m— fc 

is extreme in u m+1 ~ fc . Therefore, due to Proposition 4.7 the following matrix 

(Ak Qkxm-k ©fcxl 

Qm—kxk Qm—kxm-k (2)771— fcxl 
©ixfc (l)lxm-fc 1 
©fcxm+l-fc 

is extreme in U m+1 , and it contains the matrix A m as a principal submatrix. 

Step II. Now let us prove that U m ^ U m+1 . Let A rn e U m be any matrix. Then according 
to the Krein-Milman theorem we have 
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n . . 

where Aj > 0, Aj = 1 and Am G ExtrU m . Due to I-Step, for every i — l,m there exists a 

i=l 

matrix G ExtrU m+1 



4 (l) = 

m+l I Aki,m-h 



^ i 



'2 



containing the matrix as a principal submatrix. Here A 1 **'™ fc< = ( 0i X fc , (^), . ) , 
Then, the following matrix 



^m+l — AjA 



n 

(0 

m+lJ 

i=l 



belongs to U m+1 and it contains the matrix A m as a principal submatrix. This completes the 
proof. □ 

Corollary 4.9. We have the following inclusions: 

(i) tfcU^U^UjH-^ U"" 1 C U m _i U m C U m ; 

(ii) ExtrU 1 C ExtrUi ExtrU™" 1 C ExtrU m _! ^ ExtrU m C ExtrU m . 

5. Algebraic structure of the sets U m and U m 

In this section, we are going to study an algebraic structure of the sets U m and U m . 
Let us consider the following matrix equation 

5if ' \„, (10) 



where A m is a given symmetric matrix and X m is unknown matrix, X l m is the transpose of X m . 

In this section, we are going to solve the following problem: find necessary and sufficient con- 
ditions for A rn in which the matrix equation (10) has a solution in the class of all (sub) stochastic 
matrices. 

We will use the following result which has been proved in [4]. 
Proposition 5.1. [4] Let A rn be a extreme matrix in XJ m . If A rn has no any saturated principal 
submatrices of order m — 1 then A m is a stochastic matrix. 

Theorem 5.2. Let A m be a symmetric matrix with nonnegative entries. For solvability of 
equation 

Xm+ 2 XL =A m , (11) 

in the class of stochastic matrices it is necessary and sufficient to be A m G U m . 

Proof. Necessity. Let a stochastic matrix X m be a solution of (11). We want to show that 
A m G U m . Indeed, one can see that A l m = A m and 

Yl a *i = 2 ( ^ Xij + S x ^ ) 

i,j€a \i,jea i,jea / 

i,j£at i&ct jda ida 
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for any a C I m . Moreover, if a = I m then we have 

in | / in 



2 , 



m 



= E 



x ij 




which means that A m G U m . 

Sufficiency. Let A m G U m . We must to show the existence of a stochastic matrix X m for 
which (11) is satisfied. 

First, assume that A m G ExtrU m . In this case, we use induction with respect to the order 
of A rn . Elementary calculations show that the assertion of the theorem is true for m — 2. We 
assume that the assertion of the theorem is true for all m < k — 1 and we prove it for m = k. 

If Ak has no any saturated principal submatrices of order k — 1, then according to Proposition 
5.1, Ak is a stochastic matrix. In this case, as a solution of equation (11) we can take Ak itself. 

Let us assume Ak has a saturated principal sub matrix of order k — 1. We denote it by A a , 
where \a\ — k — 1. Since A k G ExtrU fc , due to Proposition 2.15, we have A k G ExtrUfc. Since 
A a is a saturated principal submatrix of Ak and Ak G ExtrUfc, according to Theorem 2.14 (ii), 
we get A a G ExtrUui. From \a\ = k — 1, due to the assumption of induction, for the matrix 
A a there exists a solution of equation (11) in the class of stochastic matrices. We denote this 
solution by X' a = (a^)jj eQ . Let {i } = I m \ct. We know that A k and A a are saturated matrices, 
then one can get 



+ 2 E 



Xij 



> w - = l. (12) 

jea 

By A k G ExtrUfc, according to Theorem 2.14 (ii), equality (12) yields the following possible 

two cases 

Case I: Oj j = 1 and a^j = aji = for all j G a; 

CASE II: a iojo = a joio = \ for some j ^ i and a ioj = a jio = for all j G I m \ {j }- 
In Case I, we define a solution X m = (%ij)™j =1 of equation (11), corresponding to the matrix 
A m , as follows 

% = i , j G a 

j = iq, i G a ' 

1 i = io, j = io 

where X' a = (x^)jj ea is a solution of equation (11), corresponding to the matrix A a . One can 
easily check that X m is a stochastic matrix. 

In Case II, let us define a solution X m = (xij)^ =1 of equation (11), corresponding to the 
matrix A m , as follows 

i,j G a 

% = i , j G I m \ {j } 
j = i , i G I m 
i = io, J = Jo 

where X' a = (x^)ij ea is a solution of equation (11), corresponding to the matrix A a . One can 
easily check that X m is a stochastic matrix. 

So, for extreme matrices of the set U m , the assertion of the theorem has been proved. 

Now, we are going to prove it, for any elements of the set U m . 



Xij 
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Let A m G U m . According to Krein-Milman's theorem, A m can be represented as the convex 
combination of extreme matrices of U m , i.e. 

s 

A m = J2^, (13) 



i=i 



where, < A; < 1, £ \ = 1, and ^™ e ExtrU" 1 for all i = l,s. 

i=i 

By Im, we denote solutions of equation (11), corresponding to the extreme matrices Am of 
U m , where i — 1, s. 

We define a matrix X m as follows 

s 

X m = 5>X«. (14) 
i=i 

Since every matrix Xm is stochastic, the matrix X m defined by (14) is a solution of equation 
(11) in the class of stochastic matrices, corresponding to A m . □ 

By means of Theorem 5.2 we are going to generalize this result for substochastic matrices. 

Theorem 5.3. Let A m be a symmetric matrix with nonnegative entries. For solvability of 
equation 

m \ m = A m , (15) 
in the class of substochastic matrices it is necessary and sufficient to be A m G U m . 

Proof. Necessity. Let a substochastic matrix X m be a solution of equation (15). We want to 
show that A m G U m . Indeed, one can see that A m = A m and 



°u — 2 ( 5Z 

ijeot \i,j€a 



for any ct C J m , this means that A m G U m . 

Sufficiency. Let A m G U m . We must to show the existence of a substochastic matrix X m for 
which (15) is satisfied. 

As above proved theorem, we shall prove the assertion, for extreme matrices of U m . Then, 
we prove it for any elements of U m . 

Let A m G ExtrU m . If A m is a saturated matrix, then due to Proposition 2.15, A m G ExtrU" 1 . 
According to Theorem 5.2, there exists a solution X m of equation (15), corresponding to A m , 
in the class of stochastic matrices. Since every stochastic matrix is substochastic, the matrix 
X m is a solution of equation (15) in the class of substochastic matrices. 

If A m is not a saturated matrix, then due to Corollary 4.5, there exist two index sets a 
and a' with \a'\ = K and a = I m \ a' such that A a is a saturated principal submatrix of 
A rn containing all nonzero entries of A m , and for any % G a' every i th row and i th column of 

m 

A rn consists zeroes, where K = m — Yl a ij- Since A a is a saturated there exists a solution 

m=i 

X' a = (x' i: j)ij ea of equation (15) in the class of substochastic matrices, corresponding to A a . 
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Now, using the matrix X' a , we construct a substochastic matrix X m = (xj J )™ =1 as follows 

{x'ij i,jea 
i E a', j E I m , 
j ea', ie l m 

which is a solution of equation (15) corresponding to the matrix A m . 

Hence, for extreme matrices of the set U m , the assertion of the theorem has been proved. 

For any elements of the set U m the proof can be proceeded by the same argument as in the 
proof of Theorem 5.2. □ 
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